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ABSTRACT 

Context. Large spectroscopic surveys have enabled in the recent years the computation of three-dimensional interstellar extinction 
maps thanks to accurate stellar atmospheric parameters and line-of-sight distances. Such maps are complementary to 3D maps ex¬ 
tracted from photometry, allowing a more thorough study of the dust properties. 

Aims. Our goal is to use the high-resolution spectroscopic survey Gaia-ESO in order to obtain with a good distance resolution the 
interstellar extinction and its dependency as a function of the environment and the Galactocentric position. 

Methods. We use the stellar atmospheric parameters of more than 5000 stars, obtained from the Gaia-ESO survey second internal 
data release, and combine them with optical (SDSS) and near-infrared (VISTA) photometry as well as different sets of theoretical 
stellar isochrones, in order to calculate line-of-sight extinction and distances. The extinction coefficients are then compared with the 
literature to discuss their dependancy on the stellar parameters and position in the Galaxy. 

Results. Within the errors of our method, our work does not show that there is any dependence of the interstellar extinction coefficient 
on the atmospheric parameters of the stars. We do not find any evidence of the variation of E(J - H)IE{J - K) with the angle from 
the Galactic centre nor with Galactocentric distance. This suggests that we are dealing with a uniform extinction law in the SDSS 
ugriz bands and the near-IR JHKs bands. Therefore, extinction maps using mean colour-excesses and assuming a constant extinction 
coefficient can be used without introducing any systematic errors. 

Key words. Galaxy: structure, stellar content - ISM: dust, extinction 


1. Introduction 


Understanding the dust spatial distribution in the Milky is a cru¬ 
cial part of Galactic archeology. Indeed, it can reveal important 
features of Galaxy evolution, such as the location and intensity 
of past star formation episodes ( Boulang^|2QQ7| ). Among the 
literature, the most co mmonly used full-sky dust map is that of 
[Schlegel et al.|(|1998|) obtained wit h C QBE/DIRBE data, later 
improv ed by |Schlafly et al.| ( |201Q| ) and |Schlafly & Finkbeiner| 
( |201l| ) by adding correction terms mainly due to the adopted 
reddening law. Other 2-D extinction maps result from specific 
stellar populations. For example, red clump stars are considered 
to be an ideal tracer for extinction as their mean intrinsic colour 
varies only slightly with metallicity therefore making them a re- 


* Based on observations collected with the FLAMES spectrograph 
at the VLT/UT2 telescope (Paranal Observatory, ESO, Chile), for the 
Gaia-ESO Large Public Survey, programme 188.B-300 


liable tracer of dust extinction (see for example [Gonzalez et al. 
2011[|201^|Nataf et al.|2013] ). However, such studies require a 
sufficient number of red clump stars in order to have a good spa¬ 
tial coverage and hence are limited to regions with high stellar 
densi ty (close to the plane or towards the Bulge). On the other 
hand, Nidever et al. ( |2012| ) mapped the extinction with the so- 
called Rayleigh Jeans Colour Excess method (RJCE) which is 
based on a combination of near and mid-infrared photometry 
(e.g. H and [4.5]). RJCE determines 2D star-by-star reddening 
at high-resolution (2x2') allowing to penetrate the heavily ob¬ 
scured Galactic mid-plane. 

The recent growth of large surveys together with the in¬ 
creasing volume of data has pushed forward the extinction map¬ 
ping and allowed to trace its distribution in three dimensions. 
The first of these 3D extinction models was constructed by 
Drimmel et al.| ( [200^ , by fitti ng the far and near-IR data from 
the COBE/DIRBE instrument. [Marshall et al.j ( [2006| ) used the 
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2MASS colours to gether with the ste llar population synthesis 
model of Besan 9 on ( Robin et al.|2003| ) to trace extinction in 3D 


for the Galactic Bulge region (|/| < 90, \b\ < 10) with a spatial 
resolution of 15'. 

A number of photometric techniques have been used to study 


the dust in 3D ([Green et al.[[2014|Schlafly et al.[2014[ 

Bailer- 

Jones 

2011 [[Hanson & Bailer- Jones|2014| Sale|2012[[Sa 

e et al. 

20141 

Lallement et al.|2014[).[Chen et al. 

([2014a[) traced the stel- 

lar locus similar as done by [Berry et al. ( 

2012J method by com- 


bining optical, 2MASS and WISE photometry and defining the 
reference stellar locus as well as fixing the extinction law. Their 
map covers about 6000. sq. degree with a small overlap to some 
of our GES fields (see Sect. 5) 

Until now, most of the interstellar extinction studies were 
done mainly by photometric techniques. The already available or 
upcoming large spectroscopic surveys (RAVE, APOGEE, Gaia- 
ESO , GALAH, 4MOST, etc..) provide an important quantity 
of accurate stellar parameters of different stellar populations. It 
is therefore possible to probe the 3D distribution of interstel¬ 
lar dust using the expected colours of the targets and comparing 
them with the observed ones. The Gaia space mission of ESA 
will provide in the following years two-dimensional maps for 
most of the Galaxy but also individual extinction estimates to¬ 


gether with measured parallaxes (see Bailer-Jones et al. 2013 


for more details). In this paper we will use stellar properties de¬ 
rived from the Gaia-ESO survey (GES, [Giimore et al.|2Q12| ) to 
probe the interstellar extinction in three dimensions as done by 
[Schultheis et al.| ( |2014b| ) using APOGEE data. They compared 
3D extinctio n m odels in the Galactic Bulg e region with jMarshallj 
et al. ( |2006| ) and Schultheis et al.| ( [^14a ) and found a steep rise 


in extinction in the first few kpc and a flattening of the extinc¬ 
tion atabout£^cfrom_^^ Sun. While [Schultheis et 'ST] ( [20141^ 
and Wang & Jiang [ ( [2014[ ) probed the interstellar dust proper¬ 
ties with APOGEE in the Galactic plane (|z| < 1 kpc), the GES 
fields are located at much higher Galactic latitudes (|b| > 20®). 
The GES data are thus complementary to APOGEE allowing to 
trace the dust extinction at higher Galactic height |z| and com¬ 
pare them with available 3D dust models. In Sect. 2, we describe 
the sample of stars that has been used, in Sect 3 we present the 
method employed in the determination of the extinction and the 
distances. In Sections 4 and 5, we compare and discuss the ex¬ 
tinctions to existing 2D and 3D maps, and we finish in Sect. 6 
with the discussion about the universality of the extinction law. 


2. The sample 

The Gaia-ESO survey (GES) is a public spectroscopic survey 
targeting ~ 10^ stars, covering all the major components of 
the Milky Way, from the halo to star forming regions, with 
the purpose of characterising the chemistry and kinematics of 
these populations. It uses the EL AMES multi-object spectro¬ 
graph on the VET UT2 telescope to obtain high quality, uni¬ 
formly calibrated spectra. The GES processing flow goes from 
target selection, through data reduction, spectrum analysis, as- 
trophysical parameter determination, calibration and homogeni¬ 
sation. A detailed description of the data processing cascade 
and general characterisation of the data set can be found in 
Gilmore et al. (2015, in prep.). In this paper, we analyse the 
second data release (DR2) GES results for ~ 10 000 Milky 
Way stars observed with the high-resolution gratings HRIO cen¬ 
tred at 5488 A (R ~ 19800) and HR21 centred at 8757 A (R ~ 
16200) of the GIRAEEE spectrograph. All the targets were se¬ 
lected from VISTA photometry, with colour cuts in the range 
0.2 < (J-K) < 0.8, and magnitude cuts between 12.5 < J < 17.5 


(c.f. Gilmore et al.[[2012[ ). Additional SDSS photometry is also 
available which we will use in this work. Considerable effort 
has been invested in the determination of the stellar parameters. 
The stellar parameters have been derived from three different 
methods, MATISS E ( [Recio-Blanco et al.[[2006[), SME ([Valenti[ 
[& Piskunov|1996[ ) and EERRE ( [Allende Prieto et al.[2006[ ). This 
ensures a reliable determination of the derived stellar parame¬ 
ters, a crucial step to get accurate stellar abundances and line- 
of-sight distances. The homogenisation of the results from the 
three nodes which was verified during the GES parameters vali¬ 
dation process leads to the so-called “Recommended stellar pa¬ 
rameters”. Eor our analysis we used those parameters, i.e. effec¬ 
tive temperature (Tgff), surface gravities (log g), global metallic- 
ities ([M/H]) and a-elements ([tr/Ee]). The relative error dis¬ 
tributions peak at 70 K for Tgff, 0.10 dex for logg, 0.08 dex for 
[M/H] and 0.03 dex for [tr/Ee]. More details ab out the related 
GES parameterisation pipeline can be found in [Recio-Blanco 
[erar] ( [20T4i| . 



Fig. 1. Eield centres of the GES-DR2 fields in Galactic coordi¬ 
nates 


To achieve the scope of this paper, we selected targets 
with available photometry from the VISTA Variables in the 
Via Lactea (VVV) survey ( [Minniti et al.[[2010] ) or 2MASS, to¬ 
gether with reliable atmospheric parameters obtained from spec¬ 
troscopy. This implied that we selected only the targets with 
spectra having a signal-to-noise S/N > 10, errors in radial ve¬ 
locity (t{RV) < 1.5km s“\ and low internal dispersion on Teff 
and logg between the three different nodes (see Sect. [^. The 
first two criteria ensure that the random errors of the algorithms 
are minimised, and the (crx.ff, criogg) criterium ensures a minimi¬ 
sation of the internal errors of the GES homogenisation pipeline. 
Our sample consists thus of 5603 stars with only Near-IR pho¬ 
tometry, while 1106 stars have ugrizJHK photometry. 

Eigure shows the field centres of the GES-DR2 fields in 
Galactic coordinates. 
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3. Distance and extinction determination 


The routine that has been used to determine the absolute magni¬ 
tudes of the stars in different photometric bands (and from that 
the extinctions and the distances), is based on the one described 
Kordopatis et al.](|2011b|), and already success f ully applied in 


Gazzano et al.| ( |2013| ); [Kordopatis et ^ ( |2013| ); [Recio-Blanco 
et al.| ( |2014b| ). Briefly, the method projects the measured atmo¬ 
spheric parameters (0 = Tgff, logg, [M/H]) and colours on a 
given set of theoretical isochrones. The set of isochrones used is 
defined for a given age a and a range of iron abundances [Fe/H] 
within the [M/H] error bars. For this paper, we assume that the 
[Fe/H] values of the isochrones are similar enough to [M/H], to 
use the approximation [Fe/H]=[M/H] while performing the pro¬ 
jection on the isochrones. According to Recio-Blanco et al. (in 
prep), this is true for most of the stars in the GES iDR2. 

The probability density function for a given star, 
W(a, m, [Fe/H]), is defined as: 


W(a, m, [Fe/H]) = dm • exp 




(Oj - 0i? 


( 1 ) 


where 6i is the theoretical Tgff, logg or [Fe/H], Of and cr§. 
are the measured parameters and their respective errors, and dm 
is the mass step between two points of the same isochrone, in¬ 
troduced in order to impose a uniform prior on the stellar mass. 
We note that Zwitter et al.| ( |2010 ) have shown that having such 
a fiat prior on mass does not affect significantly the final de¬ 
rived distances, when compared to the use of a combination of 
a more realistic mass function (e.g. |Chabrier|20Q3|) with a lumi¬ 
nosity prior on the surveyed stars (see, |Zwitter et al.|2010 their 
Sect. 2.2, for more details). 

The expected value of the absolute magnitude Mj in a given 
photometric band r of a given star is then obtained by computing 
the weighted mean: 


Mr 


Za,m,[Fe/H] Wju^m, [Fe/H]) • M^((2, m, [Fe/H]) 
Za,m,[Fe/H] [Fc/H]) 


( 2 ) 


where Za,m,[Fe/H] is the triple sum over the ages, masses and iron 
abundances. The associated variance of the expected absolute 
magnitude Mj is obtained by: 


(T\Mr) = 


Za,m,[Fe/H] Wia, m, [Fc/H]) • [Mr - Mria, m, [Fe/H])]2 


Za,m,[Fe/H] W(a, m, [Fe/H]) 


(3) 


We have used two s ets of iso chrones: the Yonsei-Ya le ones 
(Demarque et al. 2004| ) wit h the|Lejeune et al.|(|1998|) colou r 
tables, and the Padova ones ( [Marigo|1998[|Bressaii et ai.|2012| ). 
The Yonsei-Yale (YY) isochrones have been computed using the 
provided interpolation code of YY, from which we generated a 
set of isochrones with a constant step in age of 1 Gyr, starting 
from 1 Gyr to 14 Gyr, therefore resulting to fiat prior on the age 
of the stars. As far as the metallicities are concerned, they are 
within a range of -3 < [Fe/H] < 0.8 dex, constantly spaced 
by 0.1 dex. The ^-enhancements of the isochrones have been 
selected in the following way: 


• [Fe/H] > Odex, then [tr/Fe] = 0.0 dex 

• -0.3 < [Fe/H] < -0.1 dex, then [tr/Fe] = -rO.l dex 

• -0.6 < [Fe/H] < -0.4 dex, then [tr/Fe] = -1-0.2 dex 

• -0.9 < [Fe/H] < -0.7 dex, then [tr/Fe] = -rO.3 dex 

• [Fe/H] < -1 dex, then [tr/Fe] = -r0.4dex 


Table 1. Median errors in distance when adding an offset to the 
derived atmospheric parameters for giants (log g < 3) and dwarfs 
and subgiants (logg > 3). The “-r” sign specifies adding the 
offset, the sign subtracting the offset 



Teff + 100 K 

log g + 0.2 dex 

[M/H] + 0.1 dex 


% 

% 

% 

Giants (+) 

1.1 

19.7 

1.8 

Giants (-) 

5.1 

34.3 

1.3 

Dwarfs (+) 

5 

2.9 

7.6 

Dwarfs (-) 

11 

8.1 

9.7 


According to [Carpenter et al. (2001), the (JK)eso provided 
by the Yonsei-Yale isochrones match very well the (JKs) 2 mass 
and (JKs)vhs ^ so no colour transformation is needed when ma¬ 
nipulating the magnitudes from the different photometric sys¬ 
tems. 




Fig. 2. Left: Density plot of the difference between Teff and the 
corresponding temperature from the isochrones and the disper¬ 
sion of Teff between the three nodes (see text). The black points 
indicate the mean (in 100 K steps) and the error bars the standard 
deviations. Right: The same but for logg. The black points in¬ 
dicate the mean (in 0.2 dex steps) and the error-bars the standard 
deviations. 


As far as the Padova isochrones are concerned, they have 
been downloaded using the online interpolation interfac^ which 
allows us to select the output photometric system (2MASS, 
SDSS, VISTA). The considered metallicity range is smaller than 
the one of YY (from -2.2 dex to -1-0.2 dex in steps of 0.1 dex 
with solar-scale a-abundances), computed with steps in age of 
0.5 Gyr. 

Once the absolute magnitudes are computed, the colours are 
derived and the colour excess are deduced by subtracting the 
theoretical colour from the observed one in the five SDSS and 
3 VISTA filters. The colour excess derived using this method 
is hereafter referred to as with A = u,g, r, /, z, /, H, Kg. 

About 5% of our stars show negative extinction, most of which 
are fainter than Ks > 14.5. We omitted those from our analysis. 
The distances were calculated using the usual relation : 


logio d = 


Ks-Mk,-A[Ks) + 5 


(4) 


with d expressed in pc, and adopting A(Ks) = 0.528 x E(J - 
Kg) (jNishiyama et al.||2009 ) similar to what has been done in 
jS chultheis et al. j ( |20 14b j ) . The errors in the derived distances in¬ 
clude errors in Tgff, logg and [M/H] and errors on the appar¬ 
ent magnitude J and Kg as well as the extinction Ak, . For more 


^ http://stev.oapd.inaf.it/cgi-bin/cmd 
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Age= 5 

[M/H] = -0.2 ^ 
Age= 8 ^ 

[M/H] = -0.7 jf/ 
Age=10 

[M/H] = -1.0 

Age=12 JffY 






6000 5000 4000 3000 6000 5000 4000 3000 6000 5000 4000 3000 6000 5000 4000 3000 

Teff 


0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 0.2 0.4 0.6 0.8 1.0 

(J-K) (J-K) (J-K) (J-K) 



Fig. 3. Left: Isochrones in the Teff vs. log^ space for different combinations of [M/H] and ages. The red line shows the Padova 
isochrones. The dashed black line indicates the Yonsei-Yale isochrones without a-elements while the plain black line is for YY 
models with alpha enhancement. Right: Similar as in the left panel but in the J-K vs. K colour-magnitude diagram 


details, we refer to [Kordopatis et al. ( 2011a| ). Other than the in¬ 
ternal errors in the stellar parameters, the absolute calibration in 
Teff, log g and [Fe/H] is also crucial for the errors in our derived 
distances. 



E(B-V)Yaie distance(Padova) 


Fig. 4. Left: Density distribution of the difference in E(B-V) de¬ 
rived from the Padova and the Yonsei-Yale (YY) isochrones as 
a function of E(B-V) from the YY isochrones. Right: Ratio of 
distances derived from the Padova isochrones to the YY stellar 
library, as a function of distance. The dashed horizontal lines 
indicate +20% difference. 



-0.4 -0.2 0 0.2 0.4 

AE(B-V) 


Fig. 5. Histogram of the differences between our derived E(B- 
V) and the Schlegel SFD98 value for high galactic latitude stars 
(|b| > 10®).The black line are the derived extinction values 
using the Padova isochrones, the red line those from the Yale 
isochrones. The mean value and the r.m.s scatter is indicated on 
the top left corner. 


Table[T] shows the errors in the derived distances (in %) if one 
assumes an offset in Teff ± 100 K, logg ± 0.2 dex and [M/H] ± 



0.0 0.2 0.4 0.6 0.8 1.0 1.2 

EBV(Padova) 


Fig. 6. Density plot of E(B-V) against E(B-V) from SFD98. 



Fig. 7. Difference of E(B-V) to E(B-V) from SDF98 vs. Galactic 
latitude 


0.1 dex. Clearly seen is the large impact of logg for giants with 
up to ~ 34% error in the derived distances while temperature 
and metallicity offsets play only a minor role. For dwarf stars 
and subgiants, the effect of the surface gravity is the smallest 
(< 10%) due to the fact the main-sequences roughly overlap at 
all ages. On the other hand, offsets in Tgff have a larger effect 
for dwarf stars, because of the overall steeper slope of the main- 
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sequences in the / vs (/ - K) plane. We note, however, that the 
effect of TefF offsets rarely exceeds 10%. 

Figurej^shows the difference between the GES Tgff and log g 
and the projected value on the Padova isochrones as a function of 
(TT,ff and (Tiogg (the internal GES parameter dispersions resulting 
from the three individual nodes, see Sect. 2). It shows that when 
the internal GES parameter dispersions are large, then the ab¬ 
solute differences between the projected Teff and log g and the 
recommended GES ones also increase. In general, this indicates 
that for those few stars, the recommended parameters do not lie 
on top of theoretical isochrones and that offsets are performed 
during the projection. Given the trends illustrated in Fig.[^ we 
will use for our analysis only the targets with (Tt^^ < 300 K and 
(Tiogg < 0-3 dex, resulting to a total of 5603 stars. 


3.1. Effect of the stellar isochrone libraries 


The choice of the stellar atmosphere models used to produce 
the isochrones affect the derived extinctions and distances. In 
what follows, we compare the extinction and d istance values de- 
rived from the same pipeline (see Sect. and |Kordopatis et al.| 
201 lb| ) using alternatively the Padova an d Yonsei-Yale models. 
We recall that [Schultheis et al.| ( |2014bj ) have shown that there 
are only small differences in the derived extinction and distances 
between the P adova isochrones and the corresponding Basel3.1 
model library ( jLejeune et al.|1997] ) for K/M giants observed by 
APOGEE. Here we want to demonstrate the importance of the 
chosen model library for our GES sample covering a much larger 
Teff and log g range compared to APOGEE. 

Figure [^compares the Padova and Yonsei-Yale libraries for 
different combinations of ages and metallicities. One can see 
that the differences between the two libraries increase with de¬ 
creasing metallicity. These differences concern the positions of 
the turn-off and the giant branch, in the sense that YY has an 
offset towards cooler temperatures, i.e. towards redder colours. 
How ever, the known age-colour (or age-met allicity) degeneracy 
(e.g. [Bergemann et al.|2014| |Worthe}^1994| ) makes that the YY 
and the Padova models can partly overlap at a given metallic¬ 
ity when selecting younger YY isochrones. Since our procedure 
projects the observed Teff, logg and [M/H] on all the ages of 
a set of isochrones (see above), the differences in the derived 
magnitudes will therefore be smaller for most types of stars than 
what is suggested by a simple one-to-one comparison between 
the isochrones. However, for the stars at the boundaries of the li¬ 
braries (e.g. the hotter stars and the more metal-poor giants), the 
differences are expected to be the largest, due to the described 
offset. Finally, we also investigated the effect on the shape of the 
isochrones (and therefore on the derived absolute magnitudes) of 
atmospheric models with different ^-enhancements. The dashed 
and plain lines in Fig. obtained for the YY isochrones, in¬ 
dicate only a very small effect in the Teff vs. logg and J vs 
(J-Ks) diagram. We conclude, in agreement with other studies 
(e.g. [Breddels et al.|2010[ Zwitter et al.|2010| ), that the adopted 
(^-enhancement level of the isochrones is not affecting signifi¬ 
cantly the final distance estimations. 

Figure shows the comparison of the derived extinctions 
(left panel) and distances (right panel). It is obvious from 
this comparison that the YY E(B-V) values are systematically 
smaller than the Padova ones. In agreement with what has been 
stated in the previous paragraph, this means that the intrinsic 
colours from the YY models are redder than the Padova ones. 
The effect in E(B-V) can reach up to 0.2 mag showing that the 
choice of a certain stellar library is essential for the extinction 
determination. If one transforms this into distances (see right 


panel of Fig|^, the distances from the YY isochrones are sys¬ 
tematically larger, especially for d < 3kpc. We traced these 
differences as function of the stellar parameters, Tgff, logg and 
[M/H] and could identify that those arise mainly for dwarf stars 
(logg > 4) with 5500 < Tgff < 6500K. This is consistent with 
Fig. where the Padova isochrones predict bluer J-K colours. 
On the other hand, YY underestimate distances for cool giant 
stars with logg < 3 and 4000 < Teff < 5000K. Finally, the dif¬ 
ference between Padova and YY increases for the most metal- 
poor stars ([M/H] < -1 dex). All the above discrepancies indi¬ 
cate important differences in the stellar atmosphere models be¬ 
tween Padova and YY. For the majority of our objects, the dif¬ 
ferences between the two stellar libraries are within 20% (right 
panel of Fig. |^. In the following Section]^ we will confront the 
derived E(B-V) values with the dust map of jSchlegel et al.|19^ 


4. Comparison to the 2D extinction maps: The 
Schiegei map 


Contrary to the study of jSchultheis et al.| ( [2014b| ) with APOGEE 
where the fields were concentrated in regions of high extinction 
towards the Galactic Bulge, we analyse here w ith GES lower 
extinction fields at higher latitudes. We use the jSchlegel et al. 
{ 1998] ) dust map, hereafter referred to as SFD98, as it has the 
same sky coverage as our GES stars. We used the conversion 
E(J - Ks)/E(B - V) = 0.527 ( jRieke & Lebofsky|1985| ). 

As already mentioned, stellar libraries can give systemati¬ 
cally offset extinctions. Figureshows the difference in E(B-V) 
between the SFD98 values and the ones derived with the Padova 
isochrones (in black), and the YY isochrones (in red). We see 
clearly that the Padova isochrones match better the E(B - V)sfd98 
with a mean difference of 0.009 + 0.075 which is the typical un¬ 
certainty of our method (derived from Eq. 3). If we include low 
galactic latitude fields with \b\ <10°, the dispersion increases to 
0.18 mag where it is suspe c ted that SFD98 overestimates e xtinc- 
tion ( jSchlafly et al.||2014] ). jSchlafly & Finkbeinerj ( |2011| ) mea¬ 
sured the dust reddening using the colors of stars derived from 
stellar parameters from the SDSS. Their uncertainty is in the or¬ 
der of 30 mmag for high latitude fields and E(B - V) < 0.04. Our 
larger dispersion is due to the lower galactic latitude fields of 
the GES fields. The YY isochrones systematically overestimate 
E(B - V)sFD98 by 0.065 mag. For the remaining of our analysis, 
we decided to use the Padova isochrones, as they match better 
the SFD98 values. 

Figure shows the comparison between E(B - V)padova and 
E(B - V)sFD98- While there is no shift in the derived E(B-V), a 
large scatter is seen, especially for E(B - V)sfd98 >0.5 where 
the SFD98 values are higher than E(B - V)padova- Our results 
are in agreement with those of jSchlafly et al. ( |2014| ) where they 
compared their map based on PAN-STARRS photometry with 
the Schiegei map. They find systematically higher E(B-V) of 
SFD98 for E(B-V)sfd98 > 0-3. The Planck dust map at 353 GHz 
show a similar behaviour (see jSchlafly et al.||2014| ) indicating 
that for higher E(B-V) the far-infrared modelling of the dust 
shows some systematical offsets. A revision of these models is 
therefore needed. 


We now assess the existence of biases as a function of 
Teff, logg and [M/H] between our derived extinction and the 
SFD98 map. Figure displays a small trend of the differ¬ 
ences in E(B-V) towards higher temperatures (Teff > 6000 K) 
in the sense that SFD98 gets higher E(B-V) with respect to 
E(B - V)padova- Eor cooler temperatures (Tgff < 4500 K) SFD98 
predicts higher E(B-V) than E(B - V)padova- While for the range 
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Fig. 8. Difference in E(B-V) to SFD98 as a function of (a) Teff, (b) log g, and (c) [Fe/H] for high galactic latitude stars (|b| > 10®).The 
median value is indicated by black dots together with the standard deviation 


1.5 < log g < 4.5 no trend is visible, giants with log g < 1.5 have 
a slightly overestimated extinction compared to E(B - V)sfd98- 
Interestingly, our E(J-K) estimation for dwarf stars with log g > 

4.5 and 4000 < Teff < 6000 K is slightly offset towards higher 
values with respect to E(J - K)sfd98 • 

The derived extinction depends on the corresponding 
matched colour of the isochrone as well as of the stellar parame¬ 
ters, the surface gravity being a particularly sensitive parameter. 
Indeed, systematic offsets of 0.2 dex can significantly shift the 
extinctions and distances (see Sect. 3). This effects mainly giant 
stars while GES ha s mostly dwarf stars (see Fig. 2 of|R ecio-| 
[Blanco et al.|2014a] ). 


5. Comparison to three-dimensional maps 

In this section, we compare our 3D extinction distributions with 
the mod el of Drimmel et al. ( 2003|) and the map of Chen et al. 
( |2014bj ). The [Drimmel et aL| ( |2003j ) model is based on the far 


and near-IR data fits of the dust distribution made by [Drimmel 
[& SpergeT| ( [2001 ) from the COBE/DIRBE instrument. The spa¬ 
tial resolution of this map is approximately 21'x 21', we how¬ 
ever recall that the Drimmel et al. model does not include fea¬ 
tures related to the Galactic bar nor the nuclear disk, resulting to 
systematic overestimates of the ext inction towards the Galactic 
Bulge, (see [Schultheis et al.|2014b[ for further details). 

As far as the [Chen et al.[ ~( 2014b[ ) 3D map is concerned, 
we found only a li mited spatial overlap with six GES fields. 
Chen et ak] ( [2Q14b[ ) combined optical photometry (g, r, i) with 


2MASS (J, H K) and WISE (W l, W2) photometry and used 
the method of Berry et al. ( 2012[ ) to trace the stellar locus in a 
multi-dimensional colour space. The constructed 3D map is over 
roughly 6000 sq. degree towards the Galactic anticenter region, 
with an angular resolution varying between 3-9'. 

Figure displays for different heights above the Galactic 
plane, the Ak extinctions measured from the GES targets. The 
illustration is made in a Cartesian Galactocentric (X,Y) frame, 
with the Sun located at 8 kpc from the centre. Superimposed in 
the panel representing the closest distance from the plane, is the 
illustration of our Galaxy produced by Robert Hurt based on the 
results of the Spitzer Infrared Space telescope (R. Benjamin). 
Compared to the APOGEE targets, where the majority of them 
are located at distances larger than about 6 kpc, the GES tar¬ 


gets probe a volume much closer to the Sun with typical dis¬ 
tances d ~ 2 - 3 kpc. From Fig. one can see that GES 
also contains a line-of-sight in the direction of the Galactic bar 
(1 = 28°, b = -3°), which shows a higher concentration of dust. 
According to [Schultheis et al.[ ( [2014b[ ) this is associated to a dust 
lane in front of the bar. Clearly visible is also the increased dust 
amount associated with the Perseus spiral arm, the Sagittarius 
arm and the Scutum-Centaurus arm. Most of the low extinction 
is situated in the first few kpc around the Sun’s position. 

Figures \T0 \ and pH] show the 3D extinction for a few selected 
lines of sight and compare our results with the ones of Drim mel 
et al. and Chen et al. We present in the online table (see [A.l[ ) ad¬ 
ditional lines of sights of 3D extinction for different GES fields. 
Below we describe a few trends: 


Contrary to the APOGEE data, GES samples the first few 
kpc at higher spatial resolution in distance. In general one 
sees a steep rise in Ay (see Fig. 10 ). 

The Drimmel et al. model underestimates Ay systematically 
for high Galactic latitudes (1 > |50|). The Chen et al. map pre¬ 
dicts a steeper increase in Ay for the first kpc than Drimmel 
et al. which seem to be in better agreement with the GES 
data (see Fig. [D- 

For most of the lines of sight, we see a steep rise in Ay fol¬ 
lowed by a flattening afterwards which is qualitatively de¬ 
scribed by the Drimmel model. Due the low absolute extinc¬ 
tion values, the errors in the derived extinction and distances 
become significant for distances larger than 4 kpc. 
Puspitarini et al.[ ( [2()14[ ) studied the DIBS of 225 GES stars 


in five fields and found a good correlation between the DIB 
strength and the extinction. We have one field in common 
which is the “COROT-ANTICENTER” field (see upper left 
panel of Fig[^ located at (l,b) = (-r212.87,-2.04)^ They 
found a steep increase in extinction up to 1 kpc, a plateau be¬ 
tween 1 to 2.5 kpc and a second increase beyond 2.5 kpc (see 
their Fig. 6). We confirm the steep increase in Ay between 0 
and 1.5 kpc with a flattening starting at around 2 kpc. There 
are too few data points to see their second increase in Ay. 
The GES fields located at (l,b) = (-r38.30,-6.51)°, (l,b) = 
(-hl4.60, -r21.85)°,andthe field at (l,b) = (-r 147.13,-2.04)° 
span the full distance range and follow the Ay vs. distance 
relation predicted by Drimmel et al. 
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Fig. 9. Extinction in the (X, Y) plane for different heights above the galactic plane. On the left upper plot an illustration of our Galaxy 
produced by R. Hurt is superimposed. 


“ The GES data seem to confirm the general shape of the 3D 
extinction with a steep rise in Ay for distances up to 4 kpc 
and a flattening which occurs at shorter distance s compared 
to the values found with the APOGEE sample of [Schultheis] 
[etaL| ( |2QT4b| ). 


With the future data releases of GES in the coming years, we 
will be able to systematically trace the distance vs. Ay behaviour 
systematically allowing to compare qualitatively spectroscopi¬ 
cally derived extinction with 3D dust models. The combination 
of both GES and APOGEE data tracing the low-extinction fields 
in the Visible as well as the high-obscured fields in the Infrared 
will be clearly a goal for the future for tracing 3D extinction 


6. Interstellar extinction iaw 

Using the extinctions derived by the GES sample, we now inves¬ 
tigate the universality of the extinction law, as a function of the 
position on the sky and the stellar environment. Using APOGEE 
red clump stars, [Wang & Jian^ ( |2014| ) recently found a constant 
power law with alpha = 1.95 yielding Aj/Aks = 2.88. They 
used stars with 3500 < Tgff < 4800 K, log g < 3 and [Ee/H] > 
-1.0 dex. On the other hand, |Yuan et al.| ( [2CTT3] ) combined SDSS, 
GALEX, 2MASS and WIDE photometry to determine redden¬ 
ing coefficients from the far UV to the near and mid-IR by us¬ 
ing the SDSS spectroscopic archive. They concluded that their 
newly derived extinction coefficients differ slightly but favour 
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Fig. 10. Extinction vs distance for different lines of sight. The x-axis and y-axis give approximate the location in Galactic coordi 
gates. The Drimmel et al. (2003) model is superimposed in red 
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Fig. 11. Extinction vs distance for different lines of sight where. In red the 3D extinction model of Drimmel et al. (2003) is super¬ 
imposed while in green the 3D-model of Chen et al. (2013). 


the R(V)=3.1 Fitzpatrick reddening law ( [Fitzpatrick| 1999| ) over 
the Cardelli et al. ( Cardelli et al.|1989] ) and the O’Donnell et al. 
( Q’Donnell|1994| ) reddening laws. 


Here we perform a similar study and aim to push the 
analysis to the investigation of the systematic dependencies 
of the extinction law as a function of the stellar parameters 
(Teff,logg, [Fe/H]). Besides the near-IR photometry we use 
SDSS photometry in the filters u,g,r,i,z giving us the the pos¬ 
sibility to trace simultaneously the extinction law in the optical 
and the IR. For every band, the colour excess is simply the dif¬ 
ference between the intrinsic colour derived from the isochrone 
matching method and the observed colour. Figure shows the 
linear fit results of the different colour excess similar as done by 
[Yuan et ar] ( |2Q13| ) and |Wang & Jiang|(|2Q14|). We for ced the in¬ 
tercept to be zero as done as in |Wang & Jiang ( |2014| ). Figure 
displays our best fitting results and Table 2 gives the fitting pa¬ 
rameters as well as the comparison to |Ylian et al.| ( [2Q13| ) and 
[Fitzpatrick] ( | 1 999\ and jCardelli et al.| ( fl989| ). 


The extinction coefficients we derive agree within 10% with 
Yuan et ak] ( |2013| ) except for u-g where we find a value closer 
to the Fitzpatrick extinction law. The u-g dispersion we measure 
is higher than for other colours with some outliers not following 
the linear relation between E(u-g) and E(g-r) (see Fig. [^. The 
observed r.m.s scatter is ~ 0.1 mag in the u-band which is in 
agreement with the estimated uncertainties. 

Up to now, most of the studies have mainly been concentrat¬ 
ing on deriving interstellar extinction coefficients using a spe- 
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Fig. 12. Relation between different colour excesses. The red line 
denotes the best linear fitting. 
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Fig. 13. Extinction coefficients as a function of Teff, log g and [M/H]. 


Table 2. Derived slopes of the colour-excesses 


Colour 

This work 

Yuan et al. 

Fitzpatrick 

Cardelli 

U-g« 

0.94 +0.03 

1.08 +0.010 

0.945 

0.984 

r-E 

0.64 +0.02 

0.60 +0.010 

0.582 

0.557 

i - 

0.48 +0.01 

0.43 +0.004 

0.426 

0.496 

J-K^ 

0.41 +0.01 

0.414 +0.01 

0.411 

0.466 

J-H^ 

0.65 +0.02 

0.63 +0.01 

- 

- 


^ derived by fitting versus E(g - r) diagram 
^ derived by fitting versus E(J - K) diagram 


cific population or a mixture of stars wi th different stel lar pa¬ 
rameters. However, as first suggested by |Wilde^ ( |1963 1, other 
factors such as the stell ar abundance may alfect th e interstellar 
extinction. In particular, Grebel & Roberts] ( |1995| ) showed that 
the colour dependence of interstellar extinction is a complicated 
function of the temperature, luminosity, and metallicity of the 
stellar probe. Even the ratio of total to selective extinction de¬ 
pends on these parameters and Ry can, according to their model, 
vary more than 10% between a star having Teff = 3500 K and an¬ 
other one having Teff = 10 000 K. Eigurep^ shows the extinction 
coefficient (defined as the ratio of two colour excesses) as a func¬ 
tion of Teff and logg. It is clear from Eig. p^that our method 
does not have the required accuracy to detect variati ons at the 
level of 10% as predicted by the theoretical models of [Grebel 
[Roberts ( [1995[ ). However, within the accuracy of our procedure, 
large scale variations of the different extinction coefficients as a 
function of Tgff, log g or [M/H] are not detected. We note, on the 
other hand, that the dispersion increases for higher log g. 

Due to the limited spatial overlap between SDSS and the 
GES, only a small number of our GES sources have SDSS ugriz 
photometry. Eor the targets for which we have SDSS photometry, 
the E(g-r)/E(r-i) ratio is overestimated for cool giant stars with 
Teff < 5000 K and log g < 3.5. Within the large dispersion of our 
method we do not note any hint of a dependence of the extinc¬ 
tion coefficient to the Galactic environment, i.e. the metallicity 
of the stars. However, we lack of stars with [M/H] < -1 dex and 
small-scale variations cannot not be excluded by our method. 
Striking is the nearly flat extinction coefficient in the near-IR 
(e.g. E(J - H)/E(J - K)) plane with no variation as a function 
of the stellar parameters. Within the errors of our method, our 
work does not show that there is any dependence of the inter¬ 
stellar extinction coefficient on the atmospheric parameters of 


the stars. This suggests, that extinction maps derived from mean 
colours of stars such as the RI CE method dMajewski et al.|201 1 1 
or the colour-excesses of stars ( [Lada et al.|1994[ [Gonzalez et al.| 


2012) can be generally used, assuming a constant extinction co¬ 


efficient which does not depend on the stellar parameter. We fit¬ 


ted the E(J-H ) vs. E(J-K) diagram in the same way as [Wang ^ 
[Tiang| ( [20^ . Our value E(J - H) ^ 0.651 + 0.009 x E(J - Ks) 
is slightly higher than the one of [Wang & Jiang[ ( 2014[ ) who 
measured E(J - H) = 0.641 + 0.001 x E(J - Ks), resulting to 
a power law index of a = 2.12 and Aj/Aks = 3.15 (assum¬ 
ing Teff = 1.25 yum, 1.65 yum, 2.15 yu m for the J,H,K 5 ' bands) . Our 
power law i ndex is similar to that of Stead & Hoare (2009) with 
O' = 2.14 or [Eritz et al.[ ( [201 f] ) with a = 2.11. 

Einally, we investigated the variation of the extinction coef- 
flcient along diffe rent lines of sight. [Zasowski et al.[ ( [2009[ ) and 
Gao et al. ( 2009[ ) suggested strong variations in the extinction 


law in the mid-IR as a function of Galactic longitude or angle 
from the Galactic Center. Strong variation of the extinction law 


as a functio n of the Galactic latitu de was also found by IChen 
[et al.[ ( [2013[ ) in the Galactic Bulge. Zasowski et al.[ ( [2009| sug^ 
gested the existence of strong longitudinal variations in the in¬ 
frared extinction law where the slopes increase as the wavelength 
increase (see their Eig. 5) resulting in a steeper extinction curve 
in the outer Galaxy. [Gao et al.[ ( [2009[ ) identified small variations 
of the mid-IR extinction law with the location of the spiral arms. 
Eigure shows the extinction coefficient E(J-H)/E(J-K) as a 
function of the angle from the Galactic Center (see Eig. E). Note 
that contrary to APOGEE, GES probes different regions of the 
Galaxy and avoids especially the galactic plane where interstel¬ 
lar extinction is very high. 

Within the dispersion of our method, we do not And evi¬ 
dence for any trend of the variation of E(J-H)/E(J-K) with the 
angle from the Galactic Center nor with Galactocentric distance 
in agreement with the extinction coefficients in the SDSS bands. 
This suggests that the extinction law in the SDSS ugriz bands 
and the near-IR JHKs bands is uniform, confirming the result of 
Wang & Jiang[ ( [2014] ) obtained with APOGEE red clump stars. 


7. Conclusions 

We used data from the GES survey, together with accurate stel¬ 
lar parameters (Teff, log g, [M/H]) to trace 3D interstellar extinc¬ 
tion in intermediate and high-latitude regions of our Galaxy. We 
discuss the influence of different stellar isochrones (Yonsei-Yale 
and Padova) on the derived 3D extinction and compare our re- 
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Fig. 14. Density plot of E(J-H)/E(J-K) vs. angle from the 
Galactic Center. The dashed blanck line gives the mean of our 
derived extinction coefficient. 


suits with the SFD98 dust maps. We find on average good agree¬ 
ment with a mean difference of A£’(5 -V) = 0.009 + 0.075, with 
the dispersion getting larger when including low galactic latitude 
regions (|b| < 10°). For larger E(B - V) > 0.5 SFD98 gets higher 
extinction compared to our estimation E(B - V)padova- 

We compared our 3D interstellar dust map s with those of 
primmel et al.| ( [200^ and |Chen et al.| ( |^14b| ). The GES data 
confirm the steep rise in Ay for distances between 0 and 4 kpc 
with a flattening of the extinction at larger distances. We stud¬ 
ied the influence of the stellar parameters on the extinction co¬ 
efficients in the optical (SDSS-bands) and the near-IR (JHKs). 
We do not detect any significant dependence of the extinction 
coefficient with stellar parameters indicating that a constant ex¬ 
tinction coefficient can be assumed. Within the precision of our 
method, we do not find any evidence for the variation of the 
extinction coefficients with the angle from the Galactic centre or 
Galactocentric distance. We note, however, that our method does 
not allow to trace small-scale variations of the extinction coeffi¬ 


cient. This sug gests a uniform extinction-law, as found in [Wang 
|& Jiang| ( [2014| ). With the future data releases of GES in the com¬ 
ing years, we will be able to trace the distance vs. Ay behaviour, 
systematically allowing to compare qualitatively spectroscopi¬ 
cally derived extinction with 3D dust models 
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Fig. A.l. 3D extinctions for GES lines of sight for GES field. The red straight line gives the corresponding model of Drimmel et al. 
(2003) 
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M. Schultheis et al.: The Gaia-ESO Survey: Tracing interstellar extinction 


Fig. A. 2. 3D extinctions for GES lines of sight for GES field (continued) 
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